Large isotope effect on Tc in cuprates despite of a small electron-phonon coupling 
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We calculate the isotope coefficients a and a* for the superconducting critical temperature Tc. and 
the pseudogap temperature T* in a mean-field treatment of the t-J model including phonons. The 
pseudogap phase is identified with the d-charge-density wave (d-CDW) phase in this model. Using 
the small electron-phonon coupling constant Ad ~ 0.02 obtained previously in LDA calculations in 
YBa2Cu307, Q* is negative but negligible small whereas a increases from about 0.03 at optimal 
doping to values around 1 at small dopings in agreement with the general trend observed in many 
cuprates. Using a simple phase fluctuation model where the d-CDW has only short-range correlations 
it is shown that the large increase of a at low dopings is rather universal and does not depend on 
the existence of sharp peaks in the density of states in the pseudogap state or on specific values of 
the phonon cutoff. It rather is caused by the large depletion of spectral weight at low frequencies 
by the d-CDW and thus should also occur in other realizations of the pseudogap. 

PACS numbers: 71.10.Fd,71.10.Hf,74.72.-h,74.25.Kc 



I. INTRODUCTION 

An open and rather controversially discussed topic in 
high- Tc superconductivity is the role played by phonons 
for the low-energy electronic properties and, in particu- 
lar, the high transition temperatures T^. A direct way 
to show the involvement of the lattice in electronic prop- 
erties is the study of the isotope effect on Experi- 
mentally, the corresponding isotope coefficient a is very 
small in high- Tc cuprates near optimal doping but in- 
creases strongly in the underdoped region attaining val- 
ues being comparable or even larger than those in con- 
ventional phonon-mediated superconductors. It is of- 
ten argued that these large observed isotope effects in 
the underdoped region give direct evidence for a large 
electron-phonon coupling in the cuprates?ij^iS It even has 
been suggested that it is so large that Eliashberg theory 
breaks down and that non-adiabatic and polaronic fea- 
tures play an important role in these systems. ^i^i^ii Other 
approaches have suggested that large isotope effects may 
occur in the presence of a pseudogap.— 

Below we show that the essential features of the iso- 
tope experiments on Tc can be explained within a mean- 
field approximation of the t-J mode l^^i-'^^'^^i^^ using very 
small values for the electron-phonon coupling. The pure 
t-J model exhibits in mean-field approximation a com- 
petition of d-wave superconductivity and a c?-CDW with 
transition temperatures Tc and T*, respectively. The 
observed pseudogap can be identified with the d-CDW 
phase in this model J^ii^ Whereas many experiments sup- 
port the idea of two competing phase o^^'^^i^^'^^i^^'^°i^^ 
the nature of the additional phase remains unclear 
and many proposals besides of the c?-CDW have been 
considered,^^ Experiments suggest that its order param- 
eter has d-wave symmetry like that of the superconduct- 
ing phase. This favors an unconventional charge- or spin 



density wave state with internal d-wave symmetry rather 
than a conventional one with (anisotropic) s-wave sym- 
metry. The t-J model yields at large TV {N is the number 
of spin components) such a d-CDW but its relevance for 
the physical case TV = 2, for instance in form of a phase 
without long-range but strong d-wave short-range order, 
remains unclear j^i^liS^ 

In section II we introduce our model, its phase dia- 
gram and numerical results for the competing supercon- 
ducting and CDW order parameters, and the correspond- 
ing transition temperatures Tc and T*. We then add 
phonons assuming always that the electron-phonon in- 
teraction is very small so that they can be treated in the 
weak-coupling approximation. Explicit formulas for the 
isotope coefficients a and a* related to Tc and T* will 
be given. In section III we present numerical results for 
the doping dependence of a and a* . In section IV we ex- 
tend our treatment by including off-diagonal fluctuations 
in the d-CDW state using the method of Refs3^i22i2^. In 
this way the pseudogap phase is treated more realistically 
because the long-range order is removed and self-energy 
effects are included. Our conclusions are found in section 
V. 



II. THEORETICAL FRAMEWORK 

We consider the t-J model^^ with the Hamiltonian _ff , 



(1) 



c|^,Cio- are creation and annihilation operators, respec- 
tively, for electrons at site i and spin projection a sub- 
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ject to the condition that double occupancies of sites are 
excluded. The sums include nearest neighbor and 
next nearest neighbor {i,j)' sites on a 2D square lattice, 
the corresponding hopping elements are t and t' , respec- 
tively. Si and Ui are spin and site occupation operators, 
J the Heisenberg coupling constant, and Vc a Coulomb 
interaction between nearest neighbors. 

One way to obtain a mean-field approximation for H 
is to introduce N spin components in Eq. ([T]), scale the 
coupHng constants as t 2t/N, t' 2t'/N, J 2J/N 
etc., and to consider the large N limiti^ As a result t 
and t' are renormalized yielding the quasi-particle dis- 
persion e(k). At the same time the fermionic oper- 
ators can be treated as usual creation and annihila- 
tion operators. Explicitly, one obtains e(k) — ~2{5t 4- 
rJ){cos{kx) + cos(fcy)) — At' 5 cofi{kx) cos{ky) — fi, with 
r = X]q cos(ga;)/(e(q))- / is the Fermi function, 

S the doping away from half-filling, and fi a renormalized 
chemical potential. Here and in the following we use the 
lattice constant a of the square lattice as length unit. As 
previously discussedi^>^ the relevant order parameters 
for our mean- field treatment of the t-J model is a CDW 
order parameter 

*W = -^E^(l^-q)KT^i+QT)' (2) 
q 

with J(k) — 2J{coskx +coaky), and a superconducting 
(SC) order parameter 

^(1^) = W E("^(l^ - - ^C-lk - q))(CqtC_qi). (3) 
^ q 

Nc is the number of primitive cells, (...) denotes an 
expectation value, Q — (tt, tt) is the wave vector of 
the d-CDW, and Vc(k) = 2Vc(cosfc^ + cos ky). The 
Coulomb interaction Vc'(k) between nearest neighbors 




FIG. 1: (Color online) Zero temperature order parameters <1> 
and A and the critical temperatures T* and Tc as a function 
of doping. 



has been introduced to prevent an instability of the d- 
CDW with respect to phase separation in some region of 
phase space.— From the self-consistency condition for the 
self-energy one obtains coupled equations for $, A, the 
chemical potential and a renormalization contribution r 
to the band dispersion due to J.— Their most stable 
solutions have d-wave symmetries in the interesting dop- 
ing region, i.e., <i>(k) = $7(k) and A(k) = A7(k) with 
7(k) = (cosA:^: — cosky)/2. 

A similar mean-field approximation as above is ob- 
tained by using a slave-boson representation for H in 
Eq.©, enforcing the constraint on the average, using 
usual mean-field decouplings for the third and fourth 
term in H, and dropping the antiferromagnetic order pa- 
rameter. The above expression for e(k) as well as Eq.Q 
are in this way exactly reproduced, Eq.© with J/2 + Vc 
instead of J. 

Fig. [1] shows the doping dependence of <& and A at 
zero temperature, calculated fully self-consistently for 
t'/t = -0.35, J/t = 0.3, and Vc/J = 0.2. In the over- 
doped region (5> (5c "^0.14$ is zero. In the underdoped 
region (5 < (5c is non-zero and coexists with A. A 
shows a maximum near i5c and decays approximately lin- 
early in S towards lower and higher dopings. The two 
order parameters compete with each other which causes 
the strong decay of A with decreasing doping in the un- 
derdoped region. Also shown in Fig. [T]is T* (dash-dotted 
line) and Tc (dashed line) where $ and A, respectively, 
vanish. T* shows near 6c a reentrant behavior which 
has a simple explanation: The piece of the T* line above 
the Tc curve is unaffected by superconductivity. Discard- 
ing superconductivity the T* line would continue to the 
right decreasing slowly and reaching the x axis only at 
around 6 ~ 0.25. Taking superconductivity into account 
T* and also $ are suppressed by the presence of a finite 
A, i.e., below the Tc curve. Since A increases rapidly 
with decreasing temperature T* even bends back due to 
the strong repulsion and reaches the critical doping Sc at 
zero temperature where $ becomes nonzero. The reen- 
trant behavior thus reflects the strong competition of the 
CDW and SC order parameters. The occurrence of a 
large coexistence region of A and which extends down 
to (5 = 0, is plausible because the Fermi surface consists 
in the d-CDW state of arcs around the nodal direction^ 
which are unstable against the formation of a BCS gap 
A. 

In order to discuss the isotope effect we consider a 
phonon-induced electronic density-density coupling be- 
tween nearest neighbors and on the same atom. Approx- 
imating its frequency dependence by a rectangular form, 
as is often done in approximate solutions of the Eliash- 
berg equation, '^■^ this effective electron-electron interac- 
tion has in the d-wave channel the form, 

w(q, iujn) = -2Vnd{q, iujn)nd{-q, -iuJn), (4) 
ndiq,iuJn) = ^ 7(k)e„c^q^Ck<T. (5) 

^ k,(T=l,2 
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Similarly, we have in the isotropic s-wave channel, 
W 

w(q, iujn) = — 2""''*^' (6) 



n(q, iuJn 



1 



E 

k,cr=l,2 



(7) 



V and 14^ are electron-phonon (EP) coupling constants 
in the d and s-wave channels, respectively, 6„ the cutoff 
function Qiuo — cOn the bosonic Matsubara fre- 

quency Un = 2n7rT, and the phonon cutoff frequency. 
nd{q_,iuJn) and n(q, iw„) are electronic density operators 
with d and s-wave symmetry, respectively. Effects due 
to a small EP interaction can be taken into account in 
the curves of Fig. 1 by adding the electronic self-energy 
due to v(q, iujn) and w{q, iujn) in the form of a Fock dia- 
gram. The resulting self-consistent equations lead to an 
equation for the renormalization function Z(k, iujn) due 
to W . At T = Tc this equation can be solved directly 
yielding ZCk, iuin) = Z = 1 + As where is the prod- 
uct of W and the electronic density at the Fermi energy 
and T = Tc, i.e., it refers in general to the d-CDW state. 
A second equation is obtained which determines the SC 
order parameter A(k, ia;„), 

A(k, iun) = -4 J7(k)^ 7(k')5i2(k', zu;„0 - 



T 



4T/7(k)e„— 7(k')e„'5i2(k', zc^„0- (8) 



J is equal to J — Vc- gi2 is the 12-element of the 4x4 
matrix Green's function g. Its inverse 5"^ (k, za;„) is given 
by, 



f iujnZ — e(k) — A(k, iw„) — i$(k, ia;„) \ 

— A(k, iw„) iujnZ + eCk) i<i>(k, ia;„) 

i$(k, ia;„) icij„Z — e(k) — A(k, ia;„) 

y -i$(k,iw„) -A(k,itj„) iuj„Z + e{k) J 

(9) 

with the abbreviation k = k — Q. $(k, ia;„) is the d- 
CDW order parameter renormalized by the phonons and 
given by, 

i$(k, iujn) = -4J7(k)— Y 7(k')5i3(k', iw„') 



-H2l^7(k)e„^ Y e„'7(k')5i3(k', zc.„0, (10) 



where 513 is the 13-element of the Green's function ma- 
trix g. 

For the calculation of Tc it is sufficient to linearize the 
right-hand side of Eq. ([8]) with respect to A(k, z<jj„). Fur- 
thermore, we may neglect the phonon renormalization for 
$ in this case: Below we will be interested only in a small 
EP constant V yielding also only a small renormalization. 



Moreover, as will be shown below, this small renormal- 
ization is practically independent of the ionic mass M 
and thus may be neglected in calculating the isotope ef- 
fect on Tc- We therefore have solved Eq. pn]) using only 
the first term on the right-hand side and use the solution 
in Eq. ^ to obtain gi2. Eq. ((8|) represents an integral 
equation with two separable kernels which can easily be 
solved. Writing Eq.® as a condition for Tc we find, 



with 



{l + F^i){\ + F22)~Fl^^{), 



Fn = -2Jl dc^^^^tanhi^l 



Zoj 



'2Tc 



(11) 



(12) 



F12 — 



1 



2^c 



lUJ 

2^ 



2^Tc 



)], (13) 



and F22 = Y W JF12 ■ ip is the digamma function and 5 
denotes the imaginary part. Nd{Lo) is given by 



2lj 



^7'(k)3[G(k,c.-ir;)], (14) 



where is a positive infinitesimal quantity, 
z2-e2(k)-$2(j^) 



G(k,z) 



(z2_A2(k))(z2^Ai(k))' 



(15) 



and 



Al,: 



<k) 



2 2 



(e(k) -e(k))2 -H4$2(k). (16) 



Nd{i^) represents a weighted density of electronic states 
at Tc and is shown in Fig. [2] for two different dopings. It 
consists of a sharp peak near the energy $ due to exci- 
tations across the d-CDW gap and a structure at lower 
energies related to the van Hove singularity. It is conve- 
nient to characterize 1^ by a dimensionless EP coupling 
constant. 



Ad = VNdiO). 



(17) 



According to the above equations phonons affect Tc 
in a two-fold way, namely, via V = Xd/Nd{0) and via 
Z = 1 + Xs- Xd and As characterize the phonon-induced 
pairing interaction of d-wave and s-wave symmetry, re- 
spectively. Putting As to zero V increases Tc- To see 
this we rewrite Eq. (jlip in the form 1 + _Fii = 0. 
Fii is given by Eq. (|12p if one makes there the change 
J ^ J/(l — F12/ (1 + ^22))- This means that V increases 
J and thus increases T^. On the other hand, if we put 
V = Eq. (ini) reduces to 1 Fii = with Fn given 
by Eq. 1^ modified by J J/Z and Tc -> TcZ. Each 
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FIG. 2: Weighted density Nd{yj) 
dopings (5. 



of electronic states for two 



of these two changes diminishes T^- Thus phonons may 
lower or may increase depending which of the above 
two effects is larger. Numerical calculations indicate that 
generically the second effect dominates and that de- 
creases if one couples to phonons.^'' Most important for 
us is, however, the following observation. Our aim is 
not to determine the change in when the electron- 
phonon coupling is turned on but when the ionic mass 
M is changed. It is well known that is independent 
of the ionic mass Af , thus there will be no change in Z 
by isotope substitutions and a will always be positive. 
For small EP couplings we may even put Z = \ and keep 
only the linear term in V in the calculation of the isotope 
coefficient a — — dlnTc/dln Af. From Eq. (fTTj) one finds 
for a in this limit, 



(18) 



where the derivatives in Eq. are to be taken at the 
Tc without phonons and we also assumed lod ^ M^^'^. 



III. RESULTS FOR THE ISOTOPE 
COEFFICIENTS 

In deriving the above formulas we assumed that the 
phonon-induced interaction V affects only A but not $ 
and thus also not T*. To check this approximation we 
have calculated the isotope coefficient a* related to T* 
and defined by a* = —dlTiT*/dlTiM. The calculation 
of a* is very similar to that of a. Numerical val- 



derdoped regime are given in TABLE I for V/t = 0.04 
and uJo/t = 0.1. All values for a* are negative, i.e., a* 
shows an inverse isotope effect. However, this isotope ef- 
fect is two orders of magnitude smaller than the usual 
BCS value of 1/2 and thus tiny. Furthermore, the abso- 
lute value of a* decreases with decreasing doping quite 
in contrast to a as will be shown below. The negligi- 
ble isotope effect on T* which we found is in agreement 
with the experiment^ though there exist also data which 
have been interpreted in terms of a large isotope coeffi- 
cient a* 1^ Strictly speaking in our calculation of a the 
renormalized c?-CDW order parameter at T — Tc enters 
the density of states function Nci{uj), Eq. p^ . Since 
we have shown that T* is independent of the ionic mass 
M we may conclude that the d-CDW order parameter at 
T = Tc has also only a negligible isotope effect justifying 
the above procedure to calculate a. 

Fig. [3] shows a as a function of doping for V/t = 0.04 
and two phonon cutoffs lod corresponding to the buckling 
and half-breathing phonon modes in YBa2Cu307i^ In 
the overdoped region a is nearly independent of cod and 
S and about 0.03, i.e., very small. In the undcrdopcd 
region a monotonically increases with decreasing S and 
reaches appreciable values, for instance, 1/4 at a Tc which 
is only reduced by a factor 2 from its maximum value. 

To understand the increase of a at low dopings 
better one can rewrite Eq. (|18p approximately as 
{Nd{w)) / Nd{0){— F22) , using a low Tc approximation in 
the denominator. {Nd{uj)) is an average of Nd{uj) around 
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FIG. 3: (Color online) Isotope coefficient a as a function of 
doping for two phonon cutoffs cud- Also shown are the curves 
for $ and Tc from Fig. [T] 



TABLE I: Isotope coefficient a* for different dopings 5 



s 


0.028 


0.048 


0.064 


0.090 


0.115 


0.139 


0.164 


100 • a* 


-0.20 


-0.24 


-0.27 


-0.34 


-0.42 


-0.53 


-0.66 



ues for a* as a function of doping throughout the un- 



the phonon frequency lod over an energy interval of about 
LU]j. This interval is determined by the functions "0 in 
Eq. (fT3)) and caused by the sharp cutoff in Matsubara 
frequencies. According to Fig. [2]A^d(0) decreases rapidly 
with decreasing 6 refiecting the fact that Nd{0) is due to 
the arcs left over from the Fermi line after formation of 
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the d-CDW gap. The length of the arcs, however, de- 
creases strongly with decreasing 6. From Fig. [5] it is 
clear that for most phonon frequencies the large spectral 
weight near the d-CDW gap will substantially contribute 
to this average. As a result {Nd{o-!)) / Nd{0) 3> 1 and, 
since — i^22 is a slowly increasing function with decreas- 
ing (5, a large enhancement of a results at low dopings. 
A sharp cutoff for real frequencies, usually used in BCS- 
theory, would yield {Nd{uj)) = NdioJo)- Consequently 
a would exhibit strong resonances for ujjj $ or, more 
generally, if ujd is near well-pronounced peaks in Nd{i-j). 
Implementing the phonon cutoff in terms of Matsubara 
frequencies, as we did, corresponds to a rather soft cutoff 
in real frequencies. Such a procedure is closer to an exact 
solution of Eliashberg equations, yields a finite phonon 
contribution to Z and is also free of the above unphysical 
resonances in a if loo and the gap energy $ are of similar 
magnitude. Another advantage of our procedure is that 
a depends only weakly on lud, also for d < dc- If there 
is no pseudogap Nditu) is rather constant in the phonon 
energy region which means that the above density ratio 
{Nd{uj))/Nd{0) is one and a very small for all dopings. 

Qualitatively, the curves for a in Fig. [3] are similar 
to those in Refsi^i^ where phenomenological pseudogaps 
were used and the problem of resonances was avoided 
either by considering only the limit lojj ^ $ or by a non- 
states-conserving pseudogap. This as well as the above 
approximate expression for a in terms of the density ratio 
{Nd{u!)) / Nd{0) suggests that the above curves for a are 
rather independent of the specific features of our model 
(d-CDW with long-range order) but rather generic for 
underdoped cuprates with a pseudogap. 

Very remarkable in Fig. [3]is the fact that the tiny value 
of 0.04 for the effective EP coupling V/t is able to produce 
large values for a comparable to those seen in experiment 
in the underdoped region. The LDA yields ~ 0.022 
in YBa2Cu307,'^^ which is roughly one order of magni- 
tude smaller than X^^iS^ (For a different view on the 
magnitude of the EP coupling constants in cuprates, see 
Rel'^^). Using the LDA value Nd{0) = 1.108/el^ from 
Ref.^'', the relation Eg.pT]) and t — 0.5eV we get in the 
LDA V/t 0.04 which is the value used in our calcula- 
tion. This shows that the large experimental values for 
a in the underdoped region do not contradict, at least 
in our competing model, the small LDA values for the 
EP coupling. The case of overdoped samples is presently 
less clear because of conflicting experimental resultsjii^ 
An isotope coefficient a which is small throughout the 
overdoped region^ would agree with our Fig. [3] 



IV. EXTENSION TO FINITE CORRELATION 
LENGTHS OF THE d-CDW 

In the previous sections our employed mean field treat- 
ment yielded a d-CDW with long-range order and excita- 
tions with infinite long lifetimes. Such idealizations are 
certainly not realized in the cuprates and one may won- 



der to what degree our previous results depend on them. 
Generally speaking we do not expect drastic modifica- 
tions because the strong increase of a with decreasing 
doping was due to the rearrangement of spectral weight 
due to the d-CDW gap. This shift of spectral weight 
should not be seriously affected by fluctuations or the 
loss of long-range order as long as the correlation length 
is sufficiently large. In this section we will investigate 
this expectation on a more quantitative level. To this 
end we will employ a modelS^i^i^ which allows to treat 
exactly a certain class of semiclassical fluctuations of the 
off-diagonal order parameter. 

First we consider only the d-CDW part of g^^, i.e., 
the first and third rows and columns of Eq.Q where 
we may put Z = 1 considering again the weak-coupling 
case. Transforming the part induced by the variation of 
the order parameter into r-space we get, 

9cDwi^ ~ iV,iw„) = 

\ i7(k)$oe-*'' iujn - e(k - iV) J ' ^ ' 

with (3 = px + 4> and k = k — Q. p is a random vari- 
able with cartesian components distributed according to 
a Lorentzian. </) is a random phase which will not enter 
our final expressions and thus its distribution function 
has not to be specified. We apply now the unitary trans- 
formation 

C/i(x) = e*'3"^/^ (20) 

to Eq. (fTO|) . where 173 is a Pauli matrix. After some alge- 
bra we find, 

9cDwO^ ~ *^") = V\{^)9c\,w^ ~ iu;„)C/i(x) 
(iujn- e(k) + ei -i7(k)$o \ 
i7(k)$o i^n - e(k) + £2 / ' 

(21) 

with the abbreviations t\ = v(k)p/2 and ti = 
— v(k)p/2. Vx,y^ are given by dt{\i)/ dkx,y In de- 
riving Eq. (PT|) we also expanded e(k — iV) up to first 
order in V assuming that the changes in one-particle en- 
ergies induced by the variation of the order parameter 
vary slowly in space. To simplify the following we will 
take ti — e2 — ^ which holds exactly for the case t' — 0. 

The matrix in Eq. pT|) can be diagonalized by a second 
unitary transformation U2 yielding, 

t^l(k)5cW(k,*c^n)C/2(k) = 

iLj„-Ai(k) + e \ , , 

zc^„- A2(k)+e J ' ^^^> 

with the eigenvalues Ai,2(k) of the unperturbed d-CDW, 
given by Ea. (|16|) . Taking also superconductivity into ac- 
count we note that the Heisenberg interaction is invariant 
under the gauge transformation Ui . Applying the second 
unitary transformation U2 to the Heisenberg interaction 
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FIG. 4: Weighted density of electronic states Ndi^^) for two 
dopings S for a correlation length ^ = 250. The inset shows 
Nd{0) at low dopings. 



FIG. 5: Isotope coefficient a as a function of doping for two 
phonon cutoffs ljd calculated in the presence of short- ranged 
correlated phase ffuctuations with correlation length ^ = 250. 



and using the BCS factorization the transformed matrix 
splits into two 2x2 matrices and the resuhing gap 
equation can easily be calculated. Adding also phonons 
one finds that Eqs.(ini)-(I13l) stiU hold if the density Nd{cj) 
of Eq. (fT4| is replaced by the expression, 



2o 



k,/3 



(23) 



Finally, Eas. (fTT|) - P^ have to be averaged over the dis- 
tribution function P(p), 



(24) 



where k = 1/^ and ^ is the correlation length of the 
off-diagonal order parameter fluctuations. According to 
Eg. p^ $0 is the equilibrium d-CDW order parameter af- 
ter omitting the factor 7(k). It thus varies in r-space with 
the momentum Q because it connects electron states 
with momenta k and k. The random variable p mod- 
ulates the equilibrium total momentum Q of the d-CDW 
in an additive way and it is distributed according to a 
Lorentzian. It is sufficient to apply the necessary aver- 
age over p just to the expression of Eq. ([23|) yielding, 



1 



(c.-^/(2e)(k(k)| + |^;,(k)|))2-A2(k) 



(25) 



It is easy to see that the above expression can formally 
be obtained from Nd{u!) if one replaces there the in- 
finitesimal r] by the finite, k-dependent imaginary part 
{\v^{k)\ + \vy{k)\)/{2^). The latter quantity is in Eq. (l^ 
averaged over the Fermi line in the pure c?-CDW state. 



i.e., for (3 = 1 over the arc around the nodal line and 
for P = 2 over the remaining small piece near the 
antinodal point which however, vanishes for our two 
considered dopings. Thus the (3 = 2 contribution in 
Eq. ipS]) is negligible small. In the (3=1 contribution 
(|wj:(k)| + |t;j^(k)|)/(2^) varies only little along the arc 
so we may replace this quantity by its average over the 
arc and obtain for S = 0.085 the value 1.15vf/£, where 
vp ^ 4.6t is the square root of the Fermi surface aver- 
age of v^(k) in the normal state. This allows to describe 
phase correlations with a finite correlation length ^ as an 
inverse life time effect with energy T/t = l.lbvp/ {tS,). 

Fig. 4 shows Nd{oj) for the two dopings of Fig. 2 and 
for the case of a correlation length ^ = 250, correspond- 
ing to an inverse lifetime of about 0.02i. Though such 
a large correlation length may seem to simulate a rather 
well-ordered state most of the fine structures in Fig. 2 
are wiped out by phase fluctuations. In particular, the 
two peaks seen in Fig. 2 have merged into one broad and 
rather structureless peak. At low frequencies the changes 
introduced by phase fluctuations are rather minor. In the 
inset in Fig. 4 the static value Nd{0) is plotted as a func- 
tion of doping showing the pronounced decrease of Nd{0) 
with decreasing S similar as in the case without phase 
fluctuations in Fig. 2. This behavior for the density is 
rather robust as function of ^ as long as ^ 3> 1 holds. 

Fig. 5 shows the doping dependence of the isotope 
coefficient a using the same parameter as in Fig. 3 but 
7Vd(w) instead of iVrf(w). Comparing Figs. 3 and 5 reveals 
that this change of densities does hardly affects a so that 
corresponding curves in these two flgures are practically 
identical. This demonstrates that the steep increase of 
a with decreasing doping is not related or even caused 
by the sharp peaks present in Nd(to) or by special values 
for the phonon cutoff ujo- Instead, it is a rather univer- 
sal property caused by the large shift of spectral weight 
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towards higher frequency due to the pseudogap. If the 
correlation length ^ is decreased from large to small val- 
ues of the order of the lattice constant the depletion in 
the density Nd{0) at low energies becomes smaller and 
smaller. In accordance with the decreasing shift of spec- 
tral weight from low to high energies the isotope coeffi- 
cient a also decreases approaching a similar small value 
as in the absence of a pseudogap, i.e., in the overdoped 
region. 

Fig. 5 can qualitatively be understood in a simple 
way considering our previously approximate expression 
{Nd{uj)) / Nd{0) for a. The numerator is essentially given 
by the area under the density curve in Figs. 2 and 4. 
Its value thus is independent of the shape of the den- 
sity curve, i.e., whether it has sharp or broad peaks, as 
long as the area below the curve is constant. This area 
is practically the same in Figs. 2 and 4. On the other 
hand the oj — values of the densities both decrease 
strongly and in a similar way with decreasing S. As a re- 
sult a should be of similar magnitude in both cases and, 
in particular, show a strong increase towards low dop- 
ings in agreement with Fig. 5. The above approximate 
expression for a may also explain why the calculated val- 
ues for a of Refj^ are for our electron-phonon coupling 
V much smaller than ours: Their density of state func- 
tion at Ad = 0, plotted in their Fig. 4, is large aX uj — 
compared to the modulation due to the pseudogap which 
implies only a small redistribution of spectral weight by 
the pseudogap. 

V. CONCLUSIONS 

We have shown that a mean-field treatment of the t- 
J model which identifies the pseudogap with the gap of 



a d-CDW state is able to explain the large isotope ef- 
fect in underdoped cuprates. Interestingly, a very small 
EP coupling constant V/t 0.04 is sufficient to explain 
the experimental data. This value is very close to that 
calculated for YBa2Cu307 within the LDA. This shows 
that the large values for a of order 1 found in the un- 
derdoped region are, at least in our competing model, 
compatible with the small EP coupling constants pre- 
dicted by the LDA. The obtained huge increase of the 
isotope coefficient a with decreasing doping is rather in- 
dependent of the phonon cutoff frequency lod and the 
spectral properties of the excitations in the pseudogap 
state. The latter information is obtained by considering 
a simple phase fluctuation model where the d-CDW state 
has only short-range correlations. 



Most important for the large increase of a with de- 
creasing doping is in our calculation the large depletion 
of spectral weight at low frequencies and its shift to high 
energies by the pseudogap. Because of this we conjecture 
that our results are not specific to the employed d-CDW 
providing the pseudogap but also hold for other order pa- 
rameters such as the antiferromagnet order parameter as 
long as they lead to a strong depletion of spectral weight 
at low frequencies. Our calculation shows, in particular, 
that it is not necessary to assume a large EP interaction 
in cuprates or extrinsic effects such as pairbreaking due 
to impurities^i^ to explain the observed isotope effect 
of Tc 
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